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Abstract. The first massive stars may influence the formation of second-generation stars, in part by their metal enrichment 
of the surrounding gas. We investigate the "critical metallicity", defined as the the value, Z cr j t , at which primordial gas cools 
more efficiently by fine-structure lines of O I (63.18 jxm), Si II (34.8 Atm), Fe II (25.99 and 35.35 l±m), and C II (157.74 
flm) than by either H° or H2 line emission. We explore the time-dependent thermodynamics and fragmentation of cooling 
gas at redshifts z = 10 — 30, seeded by trace heavy elements expelled from early supernovae. Because different modes of 
nucleosynthesis (a-process, Fe-group) produce abundance ratios far from solar values, these early stellar populations are 
likely to be influenced by O, Si, and Fe cooling. Our models also include radiative coupling of the fine-structure lines and 
H2 to the cosmic microwave background (CMB), which sets a temperature floor (70-80 K at z = 25-30) that may increase 
the Jeans mass. The H2 forms from catalytic effects of electrons left over from the recombination epoch or produced during 
virialization. These electrons form the H ion (H° + e~ — ► H + 7), which in turn forms H2 through associative detachment 
(H~ + H° — » H2 + e~). In virialized halos at z = 10 — 30, the gas densities (n « 1 — 100 cm -3 ) are well below the critical 
densities, n a = 10 5 ~ 6 cm~ 3 , at which (O, Si, Fe) fine-structure lines reach LTE populations and produce their most efficient 
cooling. Thus, Z cl -[ t may initially exceed O.O1Z.0 at n ~ 1-100 cm~ 3 , and then drop to 10~ 3 ' 5 Ze, at n ~ 10 6 cm" 3 , where the 
Jeans mass may be imprinted on the stellar mass function. Primordial clouds of 10 s Mq at O.O1Z.0 and 200 K will produce 
redshifted fine-structure lines, with fluxes between 10~ 22 and 10~ 21 W m~ 2 at z ~ 4. 
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INTRODUCTION 

What is meant by the term Critical Metallicity in the 
context of first- and second-generation stars? We define 
Z cr i t as the heavy-element abundance at which metal-line 
cooling of the gas begins to dominate over cooling by H 
and He (and molecules H2 and HD). In order for the gas 
to collapse gravitationally and continue to radiate away 
the heat produced by adiabatic compression, the radiative 
cooling time, f coo i w (3nkT /2_£f), must be less than the 
gravitational collapse time, f co ii ~ (37r/32Gp) 1,/2 , where 
Jzf is the cooling rate per volume and p is the gas mass 
density. Recent studies, based on Jeans-mass arguments 
and thermodynamic histories of cloud collapse, suggest 
that the mode of star formation may change, shifting 
from higher-mass stars at low-Z (Pop III) to a more 
normal (Pop II) initial mass function (IMF) at Z > Z cr j t . 

At zero metallicity, because of the lack of CNO- 
burning and the inefficiency of the p-p chain, the first 
stars are smaller, hotter, and shorter-lived than their cur- 
rent counterparts (e.g., Tumlinson & Shull 2000). Their 
efficient production rates of ionizing radiation give them 
special importance for IGM reionization (Venkatesan, 
Tumlinson, & Shull 2003; Wyithe & Loeb 2003; Shull 
& Venkatesan 2007). The most massive stars have large 
yields of heavy elements (Heger & Woosley 2002), par- 
ticularly a-process elements (O, Si) and the iron-group. 
Thus, it is astrophysically important to understand the 
transition from first to second-generation stars when Z > 



Z CI it. However, this transition probably varies spatially 
and temporally, owing to the inhomogenous nature of 
metal production and transport into the IGM. 

In cold dark matter (CDM) cosmologies, the first 
galaxies in the universe are predicted (Ricotti, Gnedin, 
& ShuU 2002, 2007) to be 10 6 times smaller than the 
Milky Way, with characteristic masses comparable to 
mass estimates for the smaller dwarf spheroidal galax- 
ies (dSph) observed around our Galaxy and Andromeda 
(Mateo 1998; Belokurov et al. 2006). The gravitational 
potentials of these 10 6 ~ 8 M objects are so weak that 
warm and hot ionized phases of their interstellar medium 
are weakly bound. As a result, each episode of star for- 
mation may produce powerful outflows that could tem- 
porarily inhibit further star formation. 

The first subgalactic structures form from the collapse 
of rare dark matter density perturbations, with masses 
10 5_6 M Q at z ~ 30 — 40. The initial gas cooling is from 
collisionally excited H2 rotational and vibrational transi- 
tions (Lepp & Shull 1984). A minimum H2 abundance 
xh 2 ~ 10~ 4 is required to trigger star formation in a dark 
halo in less than a Hubble time. In dust-free gas, H2 for- 
mation is catalyzed by the H ion, that forms as a conse- 
quence of the shocks that partially ionize and heat the gas 
during the virialization process. At a given redshift, the 
mass of the smaller halo that can form stars is determined 
by its virial temperature and therefore by its mass. This 
analytical result has been confirmed by hydrodynamical 
cosmological simulations. 
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Abel, Bryan, & Norman (2002) carried out such nu- 
merical simulations for a selected 10 6 M halo, using 
adaptive-mesh refinement that resolves the collapse over 
a large range of scales. In this selected halo, they find 
that only one star forms, with mass 10-100 M . Bromm, 
Coppi, & Larsen (1999) found similar results with a vari- 
ety of initial conditions for the protogalaxies. These nu- 
merical results confirm longstanding theoretical sugges- 
tions that the first stars should be massive: their charac- 
teristic mass reflects the larger Jeans mass in the ineffi- 
ciently cooling metal-free gas. However, the cooling by 
trace-metal fine-structure lines depends on the gas den- 
sity (Santoro & Shull 2006, 2008) and radiative coupling 
to the CMB. Thus, the Jeans mass and critical metallicity 
are sensitive to local gas density. 

PREVIOUS CALCULATIONS OF Z crit 

The first static models of Z cr i t (Bromm & Loeb 2003) 
found Z cr ; t ~ 10~ 3 - 5 Z Q for fine-structure cooling by 
[C II] 158 jUm and [O I] 63 jUm. Santoro & Shull 
(2006) confirmed these results for C II and O I, but 
suggested that fine-structure lines of [Si II] and [Fe II] 
might also contribute, since IGM "metal pollution" from 
massive-star nucleosynthesis is weighted toward heavier 
elements. They further noted that Z cr it depends on the 
gas density, n, owing to the change in cooling rate, from 
«5f oc n 2 (at low-density) to J£ °< n (high-density), as the 
fine-structure levels reach Boltzmann (LTE) populations 
at "critical density" (n cr ). For H° excitation at 200 K, n cr 
ranges from 3000 cm" 3 ([C II]) to (1 -2) x 10 6 cm" 3 
([O I] and [Fe II]). Santoro & Shull (2006) found that 
Z cr it exceeds 0.01Z© at the low gas densities, n w 1- 
100 cm~ 3 , present in virialized halos at z > 20. 

Figure 1 shows values of Z cr ;t for C, O, Si, and Fe, 
where metallicity is labelled by a single parameter Z. In 
fact, there is no single "metallicity", since the primary 
coolants (C, O, Si, Fe) are rarely produced in solar abun- 
dance ratios. Because of the density dependence of the 
cooling, Santoro & Shull (2008) examined the thermo- 
dynamic history of cloud collapse and refined the defini- 
tion of Zcrit- In the new models they included collisional 
coupling of the gas and level populations (H2, HD, fine- 
structure lines) and radiative coupling to the cosmic mi- 
crowave background (CMB). The latter effect can be es- 
pecially important at high redshifts, z = 25 — 30, where 
the CMB temperature, 7cmb = (82 K) [(1 +z)/30], sets a 
floor on gas temperature sufficient to increase the Jeans 
mass. This, in turn, may produce more massive stars at 
high redshift (Tumlinson 2007). 



RESULTS 

The results of our time-dependent models are shown in 
Figures 2 and 3, with the astrophysical importance sum- 
marized in the captions. Figure 2 shows the thermody- 
namic (T, n) collapse history at z — 30. After an initial 
rise due to adiabatic heating, the temperature turns down- 
ward (inflection point at [Fe/H] < —2.0). Once the den- 
sity exceeds n > 10 6 cm~ 3 at Z > 10~ 3,5 Z©, T is driven 
toward 7cmb- Figure 3 illustrates the dependence of Z cr it 
on the time allowed to cool, expressed as a fraction of 
the local Hubble time. Most of the cooling time is spent 
at low densities. In order to cool in 10-15% of the gas 
must reach Z cr it ~ 0.01Z Q , until the density reaches the 
point (~ 10 6 cm~ 3 ) of most efficient cooling. 

These far-infrared fine-structure lines dominate the H2 
cooling once Z > Z cr ;t, and their detection represents a 
challenge for FIR and sub-mm astronomy. As shown by 
Santoro & Shull (2006), the strongest lines are [O I] 63 
jlim, [Si II] 34.8 £im, and [Fe II] 25.99 jim, which at 
z = 4 redshift to 316 jim, 174 jxm, and 130 jim. For 
(10 8 M Q )Mg of high-density (LTE) gas at 200 K and 
Z = (O.OlZ0)Zo.oi, the line luminosities are predicted 
to be Lii n e « ([0.8 - 2.0] x 10 41 erg s- 1 )M 8 Z .oi. The 
luminosity distance at z ~ 4 is dt ~ 10 29 cm, and the 
expected line fluxes are ~ 10~ 21 W m~ 2 , within reach of 
facilities such as ALMA, SPICA, or SAFIR. At higher 
redshifts, these lines shift into the 350 jxm (sub-mm) 
window. The fluxes will probably be lower, owing to the 
larger du lower Z, and smaller masses of cooling gas. 

OBSERVATIONS & DISCUSSION 

The "fossil record" of the first stars may be observed 
in gas-phase (IGM), as well as in metal-poor halo stars. 
Metallicities Z <~ 10~ 3 Z are measured from absorp- 
tion lines in the Lya forest at z ~ 2 — 5 (Songaila 2001; 
Schaye et al. 2003; Pettini et al. 2003; Simcoe et al. 
2004). This mean metallicity, typically inferred from 
abundances of C IV and Si IV, shows little evolution from 
z = 5 to z = 2. However, a possible redshift-dependent 
ionization correction may conspire to mask any real 
metallicity evolution. 

The origin of the metals observed in the low-density 
Lya forest at redshifts z ~ 2 — 5 is still under some de- 
bate. One view invokes a nearly uniform pre-enrichment 
of the intergalactic medium (IGM) produced by the first 
stars at high-redshift (Madau et al. 2001). The other 
view attributes the origin of the observed metal lines to 
hot, metal-enriched superbubbles located around Lyman- 
break galaxies (Adelberger et al. 2003). Given the dif- 
ficulties associated with both scenarios, it is important 
to know the amount and volume filling factor of metal- 
enriched IGM produced by primordial stars and galaxies. 
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FIGURE 1. Minimum critical metallicities, Z crlt , vs. total gas density n (Santoro & Shull 2006) for static cooling at T = 200 K 
by individual heavy elements. Curves correspond to gas enriched by C II, Si II, O I, and Fe II (see labels). Bottom envelope 
shows all four species together in solar abundance ratios. Minimum values occur at high densities (near n cr for each coolant) at 
log (Zcrit/Z©) = -3.48 (C II), -3.54 (Si II), -3.78 (O I), -3.52 (Fe II), and -4.08 (all elements). 



As modeled by many groups (e.g., Ricotti, Gnedin, 
& Shull 2002, 2007), the first sources of metals may 
come from "dwarf primordial" (dPri) galaxies, with virial 
temperatures r v ; r < 20, 000 K and circular velocities v c < 
20 km s . In contrast to more massive galaxies, the 
dark matter (DM) halos of dPri galaxies are too shallow 
to contain much photoionized gas, with temperatures 
10,000-20,000 K. During their formation, the gas is 
only heated to temperatures below 10 4 K, where it is 
unable to cool by atomic hydrogen (Lya) line emission. 
The mass of these halos is M^ m < 2 x 10 8 M Q at their 
typical redshifts of formation (z > 10). These galaxies 
rely primarily on the formation of H2 to cool and form 
stars, because metal cooling is negligible as long as the 
gas has almost primordial composition. This situation 
changes, after the metallicity rises above Z cr jt, which 
could be as large as 0.01Z Q at halo gas densities n « 
10-100 cm -3 . As the first stars form, these requirements 
no longer hold, since some gas is heated above 10,000 K 
and is polluted with heavy elements. 

At low redshifts (z < 0.4), ultraviolet spectrographs 
aboard Hubble and FUSE have made similar measure- 
ments of the gas-phase baryon content and metallicity of 
the IGM, using Lya, Lyj3, and trace metal lines of O VI, 



C III, C IV, etc. (Danforth & Shull 2005, 2007). The cur- 
rent observational sensitivity to metallicity is ~ 10~ 3 ' 5 Zo 
at high redshift (Songaila 2001; Simcoe et al. 2004) and 
1O _2 Z at low redshift (Danforth & Shull 2007; Stocke 
et al. 2007). Figure 4 shows how these lines can be used 
to derive the statistical metallicity (~ O.1Z ) from ra- 
tios of O VI to H I. However, the IGM appears to have 
widely varying metal abundances, with differences be- 
tween "filaments" and "voids" in the "Cosmic Web" of 
absorbers. If underdense regions of IGM contain some 
metals, then star formation in high redshift galaxies may 
have pre-enriched it (Madau et al. 2001). If, instead, the 
metals detected along a line of sight are associated with 
nearby bright galaxies at z ~ 2 — 5, the metal distribution 
is probably quite inhomogeneous. In this second case, 
observations are not probing a minimum floor of metal 
enrichment produced by the first galaxies. 

The degree of inhomogeneity in the IGM metal dis- 
tribution IGM is not well characterized, although this 
should change with the installation of the powerful Cos- 
mic Origins Spectrograph (COS) on the Hubble Space 
Telescope in August 2008. The spectroscopic throughput 
of COS should be over 20 x that of STIS, and the result- 
ing ultraviolet spectra will be used for a wide range of 
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FIGURE 2. Temperature-density (logT, logn) evolution of a collapsing gas cloud (Santoro & Shull 2008) starting with virial 
conditions of mass, density, and temperature at z = 30. Curves are labelled by metallicity log(Z/Z ) or [Fe/H]. The initial rise in 
temperature comes from adiabatic heating, while the subsequent decreases in T are driven by cooling from H° -excited H2 rotational 
lines and fine-structure lines of [O I], [Si II], [Fe II], [C II]. The broad temperature minima occur at n ~ 10 5,5 to 10 6 ' 5 cm~ 3 , near 
the critical densities for H° de-excitation of [O I], [Si II], and [Fe II] fine-structure levels. When level populations reach LTE, the 
volume cooling rate Jzf ' oc n rather than n 2 , and T begins to rise at n > n CI . Radiation and collisions couple the gas and fine-structure 
levels to the CMB (7cmb ~ 70-80 at z =25-30). The stellar IMF may be imprinted when T < 200 K and n « 10 6 cm" 3 . 
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FIGURE 3. Critical metallicity (Santoro & Shull 2008) for collapsing clouds, allowed to cool for various fractions (5%, 10%, 
15%) of the local Hubble time at z = 20-30. As shown by the lower two curves, the gas within virialized halos may need to reach 
metallicities as large as ~l-2% Z s in order to cool in 10-15% of the local Hubble time, f#. 



observations of the IGM and surrounding galaxies: (1) 
surveys of the baryon content and metallicity of the low- 
z IGM; (2) searches for nucleosynthetic patterns among 
heavy elements; (3) absorber-galaxy correlation studies, 
seeking connections between the IGM and the galaxies 
responsible for metal enrichment; (4) Surveys of baryons 
and metals residing in absorbers within voids. Prelim- 
inary studies of several of these issues with HST and 
FUSE spectrographs is now underway (Danforth & Shull 
2005, 2007; Stocke et al. 2006; Stocke et al. 2007). 
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FIGURE 4. The "fossil record" of the first stars may be studied in the gas phase at z < 0.4, through quasar absorption lines. 
This plot shows the "multiphase ratio" of O VI and H I absorbers, measured by Hubble Space Telescope and the Far Ultraviolet 
Spectroscopic Explorer (Danforth & Shull 2005, 2007). The ratio, N(H I)/N(0 VI) varies from high column density filaments in the 
Cosmic Web to low-column gas in voids, with Nhi < 10 14 cm -2 . Our surveys and abundance analyses show that low-z filaments 
have mean metallicities ~ O.IZ0, whereas no metals have been detected in voids, down to limits of O.O2Z0 (Stocke et al. 2007). 
The major conclusions of this study are: (1) the filaments have increased their metallicity by factors of 30-100 from z > 3 to the 
present; and (2) the IGM in voids may contain low-metallicity pockets at Z < 0.01Z Q . 



Critical Metallicities for Second-Generation Stars 



February 1, 2008 



6 



